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Ah;  TRACT 


Methods  for  specifying  pre-  and  post-compensa- 
tion to  achieve  a desired  sysiea  type  in  linear  rnulti- 
variable  servomechanisms  are  developed.  Sufficient 
conditions  are  derived  for  both  forms  of  compensa- 
tion when  the  compensator  is  restricted  to  a diago- 
nal input-output  representation,  and  a method  for 
realizing  reduced-order  pre-eonpensators  is  present- 
ed for  the  case  where  the  restriction  is  removed. 

I.  INTRODUCTION 

The  basic  theory  of  system  type  for  linear 
multivariable  (ta~ input /ra-output)  servomechanisms  has 
been  developed  in  [l-4].  There,  system  type  is  ex- 
pressed as  an  m-vector  of  integers  which  describes 
the  steady-state  tracking  properties  of  each  input- 
output  pair  of  the  closed-loop  system  as  a function 
o2  certain  characteristics  of  the  open-loop  system. 

In  order  to  exploit  this  theory,  it  is  necessary  to 
h3ve  systematic  compensation  methods  which  can  be 
used  to  change  a system  of  a given  type  to  one  of  a 
different  (higher)  type  to  satisfy  performance 
specifications.  In  a logical  sequel  to  [l],  the  pre- 
sent paper  addresses  this  issue  of  compensat ion . 

Some  votk.  along  these  lines  has  been  presented  in 
[4],  but  the  work  here  provides  a more  general  and 
comprehensive  treatment. 

In  contrast  to  the  single-input/single-output 
case  where  rendering  a system  type  L is  relatively 
straightforward  [5],  the  multivariable  case  requires 
consideration  of  both  pre-compensation  and  post-com- 
pensation. This  is  due  to  the  fact  that  rational 
transfer  function  matrices  do  not  commute  whereas 
scalar  transfer  functions  do.  A distinction  is  thus 
made  throughout  the  presentation  between  these  two 
forms  of  compensation. 

Briefly,  the  content  of  the  paper  is  the  follow- 
ing. In  Section  II,  results  are  developed  which  per- 
mit a designer  to  change  system  type  using  simple 
diagonal  transfer  function  matrix  compensators. 
Sufficient  conditions  for  attaining  any  specified 
type  with  both  pro-  and  post-compensators  are  given. 
Also  in  this  section  is  a simplified  proof  of  the 
controllability  condition  of  [b,  7]  for  a 

•Research  supported  in  part  by  the  U.S.  Air  Force 
Office  of  Scientific  Research  under  Grant  No.  AFOSR 
71-21151)  and  in  part  by  the  National  Science  Founda- 
tion under  Grunt  No.  OK  19715. 


multivariable  servomechanism  which  employs  pre- 
compensation.  Section  III  deals  with  reduced-crder 
pre-compensation.  The  main  risult  there  is  an  exist- 
ence theorem  whose  proof  is  constructive.  The  con- 
clusion is  given  in  Section  TV. 

II.  SIMPLE  COMPENSATORS 

The  first  method  of  creating  an  open-loop  trans- 
fer function  matrix  of  specified  type  is  the  simplest 
and  is  called  pcst-compensatlon.  The  reason  for  the 
name  is  obvious  from  Fig.  1. 


Fig.  1.  Post  -co«ep«nsat  ion. 

If  H (s)  is  restricted  to  be  diagonal,  the  solu- 
tion of  tfie  compensation  problem  is  characterized  as 
follows. 


Proposition  1.  Let  G(s)  be  an  m x m type 
tj  •••  transfer  function  matrix.  In  order 

that  tne  compensated  systemf  Gc(s)  * Ilc(s)  G(s)  be 

type  [ll  ll  ...  E']  where  E.  > E . for  every  i, 

1 c m 11 

1 < ( < a,  it  Is  sufficient  that 


Hc(s) 


Proof.  With  the  given  Hc(s), 

G (s)  - H (s)  G(s) 
c c 


- H (s)  H(s)  G'(s) 
c 


where,  from  the  definition  [l,  2]  of  a type 

[£,  £„  ...  I 1 transfer  function  matrix, 

1 2 m _ - 


H(s) 


ding 


1 


1 


s 


and  G1 (s)  has  no  zeros  at  the  origin.  Thus 

G (s)  “ U(s)  C'(s)  with 
c 


606 


li(n)  - li  (.:)  H(a) 


so  that  Lh»  compensated  system  is  ot  the  desired 
type.  □ 

In  gcneril*  this  configuration  requires  t he 
fewest  number  oi  integrators.  However,  using  the 
results  ot  Section  ill,  it  can  be  shown  that  there 
are  ca- * where  preco;:;»viis«it  ion  ot  lower  order  exists. 

Due  to  the  nature  of  the  outputs  of  tn  «ny  sys- 
tems, post-Co- ,»ensat icn  may  be  undesirable  or  even 
impossible  to  implement.  In  this  case  one  may  use 
pr**-corpensat  ion  as  indicated  in  Fig.  2 to  produce  a 
transfer  fuuction  matrix  of  specified  type.  The 


-CO  0 


0 0 0 


B » [B  0 0 


C . [0 


Pr*- 

Coapensat^r 
H (s) 


Fig.  2.  Pre-coaipejj«3tlon. 

structure  here  is,  of  course,  identical  to  that  of 
classical  servomechanism  theory.  The  following  pro- 
vides a sufficient  condition  for  this  configuration 
when  H (s)  is  again  restricted  to  be  diagonal. 

Proposition  2,  Let  G(s)  be  an  m x m type 

[£,  t ~ ...  i "j  transfer  function  matrix.  In  order 
12  m 

that  the  compensated  system  G (s)  ® G(s)  H (s)  be  at 

» i i ^ • c 

least  type  [£.  ...  I ] where  l.  > l.  for  at  least 

'12  m i i 

one  i,  1 < i < m,  it  is  sufficient  that 


1_  i__  1_ 

L L L 

s s s 


where  H is  m x m and  L * max  (l.  - £.). 
c ^ 1 i 

Proof . The  desired  result  is  an  immediate  con- 
sequence of  Theorem  4 of  [l]  with  G(s)  ■ G. (s)  and 
Hc(s)  - G2(s).  □ 1 

This  is  the  form  of  compensator  which  was  used 
by  Young  and  Willems  [6],  and  Porter  and  Bradshaw 
[7],  Their  approach  employed  the  state-space  de- 
scription of  the  plant  with  state  variable  feedback. 
Although  not  explicitly  so  stated,  their  solution 
was  based  on  the  fact  that  the  plant  was  at  least 
type  [0  0 ...  0]  since  their  final  condition  for 

controllability  is  equivalent  to  the  plant  having  no 
zeros  at  the  origin. 

The  theorem  given  below  is  an  extension  of  the 
basic  result  in  [6].  The  present  proof  is  new  and 
■uch  simpler  than  the  one  in  [7]. 

Consider  the  system 

x ■ Ax  + Bu(t)  + Gy^(t)  (l) 


* - [x  *2  zr^ 


where  superscript  T denotes  the  transpose.  In  these 
relations.  A,  B and  C are,  respectively,  n x n,  n x r 
and  n x n;  the  z are  E-vectors;  u(t)  is  an  accessible 
input  to  the  system;  y^ft)  is  an  m- dimensional  command 
input;  and  I and  0 are,  respectively,  the  identity 
matrix  and  null  matrices  of  suitable  order.  It  is 
shown  in  [?]  that  if  there  exists  an  r x (p  + mr) 
matrix  K such  that  A + BK  has  open  left  half-plane 
eigenvalues,  then  for 


where  t > 0 and  the  are  arbitrary  real  constants, 
the  tracking  error  e(t)  * y^ft)  " y(t)  “ yd (t)  - 
Cx(t)  has  the  property  that 

lim  e(t)  - 0 (2) 

t— » 

The  resulting  m-input/m-output  closed-loop  system 
thus  becomes  a type  [r  r . . . r]  multivariable 
servomechanism.  It  is  clear^that  the  condition  in 
Eq.  (2)  can  be  achieved  if  (A,  B)  is  controllable. 

Theorem  1.  (A,  B)  is  controllable  if  and  only 

if 

(i)  (A,  B)  is  controllable 


(it)  rank 


r A b 1 

I - n + m 

L-c  oj 


Proof ■ from  [8],  (A,  B)  is  controllable,  if  and 
only  if 

rank  [si  - A j B]  ■ n + mr  (3) 

for  all  s.  Expanding  the  matrix  In  Eq.  (3), one  has 


1* 


r 


»I  -A  0 0 ...  0 


hi  0 ...  0 0 0 


0 -I  si  ...  0 


[si  - a ' a] 


0 . . . rI 


0 . . . -I  hi 


For  s + 0,  the  la:>t  wr  rows  are  clearly  Independent, 
and  independent  ol  the  first  n rows.  Therefore,  if 
rank,  [si  - A j Fij  * n,  Eq.  (3)  would  he  satisfied 
for  s / 0.  However,  this  is  the  case  by  contl i lien 
(i).  Conversely,  if  condition  (i)  lr»  not  satisfied, 
rank  [si  - A j h]  is  strictly  les.->  than  n and  t lie 
system  is  not  controllable. 

For  s *■  0,  the  last  (n  - L)r  rows  are  independ- 
ent, and  independent  of  the  first  n -f  ui  rows.  By 
condition  (ii),  the  first  n + n rows  are  independent 
so  that  tq.  (3)  is  again  satisfied.  Conversely,  if 
condition  (il)  is  not  satisfied,  then  Eq.  (3)  will 
not  hold  for  s *•  0. 

Thus  Eq.  (3)  holds  if  and  only  if  conditions 
(1)  and  (ii)  are  tnet.^  This  is  equivalent  to  the 
controllability  of  (A,  B) . Q 

The  fact  that  Eq.  (1)  describes  a type 
[r  r . . . r]  system  can  also  be  shewn  via  the 
application  of  Theorem  3 of  [l]  and  Proposition  2 
above . 


III.  REDUCED-ORDER  PRE-COMPENSATION 

The  results  of  the  preceding  section  show  that 
unity  feedback  systems  of  a specified  type  may  be 
synthesised  using  diagonal  pre-  or  post-compensation. 
The  sufficient  condition  of  Proposition  2 gives  an 
upper  bound  on  the  order  of  the  dynamics  required  by 
a pre-compensator  to  create  the  desired  type.  For 
post-compensation.  Proposition  I showed  that,  in 
general,  a compensator  of  lower  order  could  be  used. 
In  this  section,  a technique  for  pre-compensator  de- 
sign will  be  presented  for  the  case  whete  the  re- 
striction to  diagonal  compensation  is  lifted.  The 
order  of  the  compensator  will  be  no  greater  than  the 
order  of  the  post-compensator  of  Proposition  1. 

Mathematically,  the  problem  of  pre-compensator 

design  to  create  a type  [«!.  ...  £ 1 transfer  func- 

i c n 

tlon  can  be  described  as  follows:  Given  a plant 

transfer  function  G(s),  find  a dynamic  compensator 
with  transfer  function  K(s)  so  that 

G(s)  K(s)  - H(s)  C’(s) 


H(s)  « dlag 


fj i_  i 

s <1  V * '<« 

. S 8 H 


and  G*  (h)  h.is  no  zero*  at  the  origin. 

The  technique  lor  solving  this  problem  with 
reduced-order  dvnimlen  will  he  presented  in  the  con- 
structive proof  to  theorem  2 which  is  the  main  re- 
sult of  this  section.  The  idea  of  the  theorem  is 
that  to  increase  any  of  the  by  one  in  a type 


[i‘j  l ...  transfer  (unction,  a pre-compensator 

may  be  used  which  lias  only  one  finite  pole,  i.e., 

K(h)  can  be  Found  with  order  equal  to  one.  The  proof 
of  the  theorem  requires  the  following  two  lemmas. 

Litodiu  1 . Any  column  of  an  m x m polynomial  matrix 
l>(s),  which  has  elements  which  are  riot  divisible  by 
the  linear  term  (s  - a),  can  be  brought  to  a form 
such  that  every  element  of  the  column  except  one  is 
divisible  by  (s  - u)  via  multiplication  by  a constant 
nonsingular  upper  triangular  matrix  M which  operates 
on  the  left  (i.e.,  a row  operation  matrix). 

Proof . Without  loss  of  generality,  consider  a 
single  column  of  b(s),  any  column  q,  which  has 
elements  d^(s)  In  row  i.  by  assumption  not  all  of 
the  d.(s),  i - i,  2,  ...,  m,  are  divisible  by  (s  - a). 
Therefore  there  exists  a greatest  index  k such  that 
cj^(s)  is  not  divisible  by  (s  - ex).  Polynomials  are 

elements  of  a euclidian  ring  [10]  so  d^(s)  can  be 
written  as  <^(s)  “ (s  - u)  c^(s)  4-  where  c^(s)  is 
a polynomial  in  s and  r^  is  a constant  tern*.  Also 
for  every  j < k,  d^(s)  can  be  written  as 
dj(s)  * (s  - a)  c^(s)  4-  where  each  c^(s)  is  a poly- 
nomial and  each  r.  is  a constant.  Now  if  a constant 
J 

nonsingular  upper  triangular  matrix  M is  formed  which 
operates  on  the  left  of  D(s)  so  that  -(r./r,  ) times 

J R 

row  k is  added  to  each  row  j for  j < k,  the  resulting 
polynomials  in  row  J of  column  q are 

d,(s)  * d^(s)  - (r^/rk)  dfc(s) 

- (s  - a)c^(s)  + - (r^/r^Hs  - cOcfc(s)  - 

- (s  - a)  [c j (s)  - (r7rk> 

Thus  3^(s)  is  divisible  by  (s  - a).  By  assumption, 

(s  - a)  divides  dj(s)  for  i > k (M  leaves  row  i, 

1 > k,  unchanged)  and  after  operation  on  the  left  by 
M,  (s  - a)  divides  dj(s)  for  j < k.  Thus  only  ele- 
ment d^fs)  of  column  q is  not  divisible  by  (s  - a).  □ 

Lemma  2.  Let  the  m X m polynomial  matrix  F(s) 
be  the  Smith  form  [8]  of  the  m X m polynomial  matrix 
G(s)  so  that 

F(s)  - M(s)  C(s)  N(s) 

where  M(s)  is  a unlmodular  m x m polynomial  matrix 
representing  row  operations  and  N(s)  is  a unimodular 
m x m polynomial  matrix  representing  column  opera- 
tions, and  M(s)  and  N(s)  are  not  unique.  Then  for 
every  unlmodular  row  operation  matrix  M'(s)  of  the 
type  which  adds  a polynomial  multiple  of  row  1 to  row 
J with  i > J,  there  exists  a unlmodular  column  opera- 
tion polynomial  matrix  S'(s)  such  that 

M'(s)  F(s)  N’(s)  - F(s) 

Froof . By  definition,  F(s)  is  a diagonal  poly- 
nomial matrix  with  elements  tj(s)  in  the  (1,  1)  posi- 
tion. The  Lj(s)  are  such  that  r^(s)  divides  ck<s) 

whenever  J < k.  Thus  If  M'(s)  Is  a row  operation 
polynomial  matrix  which  adds  p(s)  times  row  k to  row 
j with  | < k,  tile  result  Is 


n 


C08 


D(r)  » K P(s)  lj(s) 


M'(s)F(s) 


tj(s)  0 0 (p(s)<k(R))  0 


V8)  ° 


Since  divides  €^(s),  the  column  operation 

Datrix  which  adds  (-p(s)  r.  fs)/c^(s))  times  column  j 

to  column  V.  is  a polynomial  matrix  and  gives  the 
desired  result.  □ 

It  is  remarked  that  by  repeated  application  of 
this  lemma,  the  result  holds  ior  any  upper  triangular 
row  operation  matrix  M"(s). 

The  main  result  of  this  section  can  now  be 
proved . 

Theorem  2.  Given  an  o x c type  ft.  £„  ...  / 1 
1 2 nr 

transfer  function  matrix  G(s),  then  for  any  i, 

1 < i m,  there  exists  an  n x m transfer  function 
matrix  K(s)  of  order  equal  to  one  such  that  the 
series  r.  x m transfer  function  natrix  C(s)  K(s)  is 


type  [£.  I,  ...  t.  + 1 ...  4 ]. 

Proof . The  proof  is  constructive.  First, 
factor  G(s)  into 

C(s)  - H(s)  G'(s)  - H(s)  D_1(s)  S(s) 


[+>+>■ 


G' (s)  is  type  [0  0 . . . 0],  and  D(s)  and  N(s) 

are  polynomial  matrices  obtained  by  a coprime  factor- 
ization of  C' (s) . 

Now  the  problem  may  be  restated  as  follows: 

Find  an  o x m transfer  function  K(s)  so  that  the 
series  transfer  function  matrix  G'(s)  K(s)  is  type 
[0  0 ...  0 1 0 . . . 0]  where  the  1 is  in 

an  arbitrary  specified  position  i.  By  Theorem  1 of 
Cl] . this  is  equivalent  to  finding  K(s)  such  that  the 
coprine  factorization 


C'(s)  K(s)  - D'_1 (s)  N’(s) 


is  such  that  every  element  of  column  i of  D'(s)  is 
divisible  by  s. 

The  construction  of  K(s)  proceeds  as  follows. 

Let 

N(s)  - P(s)  N(»)  Q(s)  (5) 

where  N(s)  is  the  Smith  form  of  N(s);  and  P(s)  and 
Q(s)  are  unimodular  polynomial  matrices.  By  the 
maximal ity  of  the  there  exists  at  least  one  ele- 
ment in  every  column  of  D(s)  from  Eq.  (4)  which  is 
not  divisible  by  s.  With  P(s)  defined  by  Eq.  (5), 
find  an  upper  triangular  constant  m x ra  matrix  R so 
that  every  clement  of  column  t of  the  product 
P(s)  P(s)  is  divisible  by  s except  one.  This  is 
always  possible  by  Lemma  1 with  it  - 0.  Let  k indi- 
cate tJu  row  index  of  the  element  of  column  i which 
is  not  divisible  by  u,  and  let 


where  every  element  of  column  1 of  D(e)  in  divisible 
by  s except  the  on«*  In  row  k. 

Now  find  a uni modular  polynomial  matrix  S(s)  so 

that 

R N(s)  S (s)  « N(r.)  (7) 

Tills  is  always  possible, by  Lemma  2 and  the  remark 
that  followed  it  since*  N(s)  is  in  Smith  form.  In 
fact,  S(s)  Is  determined  uniquely  by  K. 

Using  Eq.  (6),  D * can  be  written 

b_1(s)  - frJ(s)  [k  P(s)]  (8) 

Then  with  kqs.  (5)  and  (7),  N(s)  can  be  written 

N(s)  - [R  P(s)]-1  N (s)  [Q(s)  S(s)]"1  (9) 

Now  construct  K(s)  as 

K(s)  = [Q(s)  S(s)]  K(s)  (10 


K(s)  - dlag  [1  1 ...  1 i 1 ...  l] 

with  the  1/s  factor  In  the  kth  position.  Since  both 
Q(s)  and  S(s)  are  polynomial  matrices,  K(s)  has  order 
equal  to  one. 

From  Kqs.  (8),  (9)  and  (10),  it  follows  that 
G' (s)K(s)  - D_1(s)  N(s)  [Q(s)  S(s)]  K(s) 

- D'1(s)[RP(s)][RP(s)]‘1N(s)[Q(s)S(s)]"1[Q(s)S(s)]K(8) 

- D_1(s)  N(s)  K(s>  (11) 

Since  both  K(s)  and  K(s)  are  n x m diagonal  ma- 
trices whose  elements  are  from  a commutative  ring, 
the  matrices  commute  and  Eq.  (11)  becomes 

G'(s)  K(s)  - D_1(s)  ii(s)  fi(s) 

- [K‘1(s>  D(s)]'1  N(s) 

- D_1(s)  N(s)  (12) 

By  construction  of  K(s),  K *(s)  Is  a polynomial 
matrix  which  multiplies  row  k of  D(s)  by  the  factor 
s.  Thus  every  element  of  column  1 of  D(s)  is  now 
divisible  by  s.  , „ 

It  remains  to  show  that  D(s)  and  N(s)  in  Eq.  (12) 
are  coprime.  Then  by  Theorem  1 of  [l],  the  present 
theorem  is  proved.^ 

Now  D(s)  and  N'(s)  are  coprime  if 

rank  [6 (s)  N(s)]  - m (13) 

I 

for  all  s [8] . 

Consider  first  the  case  where  s - 0.  Since  N(s) 
has  no  zeros  at  the  origin,  Eq.  (13)  Is  obviously 
true  for  s ■ 0. 

, Next  consider  the  case  where  s 4 0.  Factor 
[D(s)  i N(r)]  into  the  following: 

[D(a);N(s)]  - [if1 (s)RP(s)n(s)  j RP(s)N(s)Q(s)S(s)] 

- K*1(9)[RP(a)]l0(s)!N(s)]  “ |k'(r) 

10  Q(s)S(s)  J 
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where 


2 


0 


0 


0 Q(«)S(a) 


is*  a 2u.  x 2r»  raatrix  with  the  0;  x m block  Identity  I 
and  the  2 u x 2ai  matrix 

K'(s)  “ d lag  [1  l ...  i 1 ...  l] 

has  the^l/.s  factor  In  the  m + k position.  For 
6 i 0,  K"^(s)  and  K* (s)  are  nonsinguiar.  Therefore, 


[D(s) |N(s)] = [RP(s)]"1K(s)[D(s):  N(s)]k,-1(s) 


[0  Q(s)S(s)  J 

and  it  follows  that  rank  [o(s)'N(s)]  S rank 
[D(s)|  N 1 s ) ] = m where  the  last  equality  is  from  the 
definition  of  the  coprime  factorization  in  Eq.  (4). 

Thus  Eq.  (13)  is  satisfied  for  all  s / 0 and 
hence  for  all  s.  This  completes  the  proof.  □ 

The  proof  of  this  theorem  shows  explicitly  how 
to  find  K(s).  The  construction  nay  be  complex  since 
computation  of  the  Smith  form  for  large  systems  can 
be  difficult.  Whether  or  not  K(s)  can  be  found  as  a 
proper  rational  matrix  is  an  open  question. 

Example.  Consider  the  problem  of  determining 
the  pre-compensator  which  will  render  the  type 
[0  0 0]  transfer  function  matrix 


type  [0  0 l]. 


First  G(s)  is  factored  into  the  form 
G(s)  - D“hs)  N(s)  yielding 


b 7 (s+3) 


-2  (1-s) 


0 -7s  -9s-4 


Next,  P(s)  D(s)  Is  determined  to  be 


P(s)  D(s)  - 


- 1 7s+9  Xs+1 ) -3  (s+1) 


For  this  example.  It  is  immediately  seen  that  K is  the 
identity,  and  consequently,  that  S(s)  is  also  the 
identity.  So  K(s)  diag  [l  1 1/s].  Then,  K(s) 

is  found  to  be 


-2 


7 (s+3) 


Then  N(s)  can  be  found  as  N(s)  “ P(s)  N(s)  Q(s), 


K(s)  ■=  Q(s)  K(s)  - 0 


so  that  G(s)  F.(s)  Is 


G(s)  K(s)  =2  0 0 

6 7s+9  -7s2-9s-4 

_s+3  s+3  s(s+3) 

which. Is  of  the  desired  type. 

IV.  CONCLUSION 

As  a sequel  to  [l],  this  paper  has  presented 
methods  for  specifying  pre-  and  post-compensation  to 
achieve  a desired  system  type  for  linear  multivariable 
servomechanisms.  Sufficient  conditions  have  been  de- 
veloped for  both  forms  of  compensation  when  the  compen- 
sator is  restricted  to  be  characterized  by  a diagonal 
transfer  function  matrix.  In  the  case  of  pre-compen- 
sation, a method  of  determining  a reduced-order  com- 
pensator has  been  derived  and  illustrated. 

Further  research  Is  needed  to  extend  the  present 
results  to  permit  the  realization  of  compensators 
which  simultaneously  achieve  a desired  system  type, 
guarantee  closed-loop  stability,  and  satisfy  other 
performance  specifications  such  as  minimum  integral- 
square-error.  Some  results  in  this  direction  are  re- 
ported in  [ll]. 
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